Quantum Tera-Hertz electrodynamics in Layered Superconductors 
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In close analogy to quantum electrodynamics, we derive a quantum field theory of Josephson 
plasma waves (JPWs) in layered superconductors (LSCs), which describes two types of interacting 
JPW bosonic quanta: one massive and the other almost-massless. We also calculate the amplitude 
of their decay and scattering. We propose a mechanism of enhancement of macroscopic quantum 
tunneling (MQT) in stacks of intrinsic Josephson-j unctions (SIJJs). Due to the long-range interac- 
tions between many junctions in the LSCs, the calculated MQT escape rate T has a very nonlinear 
dependence on the number of junctions in the stack. This allows to quantitatively describe striking 
recent experiments in Bi2212 stacks. 

PACS numbers: 74.72.Hs, 74.78. Fk, 



The recent surge of interest on Stacks of Intrinsic 
Josephson Junctions (SIJJs) is partly motivated by the 
desire to develop THz devices, including emitters [llLi, 
niters Q , detectors Q , and nonlinear devices || . Macro- 
scopic quantum tunnelling (MQT) has been, until re- 
cently, considered to be negligible in high-T c supercon- 
ductors due to the d-wave symmetry of the orderpa- 
rameter. Recent unexpected experimental evidence [3,13 
of MQT in layered superconductors (LSCs) could open 
a new avenue for the applicability of SIJJs in quantum 
electronics This requires a quantum theory for SI- 
JJs capable of describing quantitatively this new stream 
of remarkable experimental data. In contrast to a single 
Josephson junction, SIJJs are strongly coupled along the 
direction perpendicular to the layers. This is because the 
thickness of these layers is of the order of a few nm, which 
is much smaller than the magnetic penetration length. 
This results in a profoundly nonlocal electrodynamics 

2] 

that strongly affects the quantum fluctuations in SIJJs. 

Using a general Lagrangian approach, we derive the 
quantum electrodynamics of JPWs, which describes two 
interacting quantum fields. We analyze the first-order 
Feynman diagrams for: (i) the decay of a quantum JPW 
propagating along the layers and, (ii) JPW- JPW scatter- 
ing. Employing the quantum statistics of these plasmons, 
we calculate the average energy of the JPWs as a func- 
tion of temperature, and find it to be much higher than 
for the same number of non-interacting junctions. Using 
this general approach, we develop a quantitative theory 
of the MQT in SIJJs. For example, we derive the MQT 
escape rate, T, which is strongly non-linear with respect 
to the number of superconducting layers, N, and changes 
to r oc N when N exceeds a certain critical value N c . Our 
results are in a good quantitative agreement with recent 
very exciting experiments . 

Quantum theory for layered superconductors. — The 
electrodynamics of SIJJs can be described by the La- 



grangian: 
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gauge-invariant interlayer phase difference, and p n = 
( s / ^ab)d x Xn — Zir^sA^f 1 /<& is the normalized supercon- 
ducting momentum in the nth layer. Here, we introduce 
the phase Xn °f the order parameter, the interlayer dis- 
tance s, the in-plane X a b and out-of-plane A c penetration 
depths, the anisotropy parameter 7 = \ c /\ a b, flux quan- 
tum $0, and vector potential A. The in-plane coordi- 
nate x is normalized by A c ; the time t is normalized by 
X/u>j, where the plasma frequency is wj; also, d x = d/dx, 
dyfn = A q6 (/„ + i - f n )/s, and ' = d/dt. We choose the z 
axis pointed along the magnetic field. Varying the action 
S = Jolt C produces the dynamical equations 



Cp n - d^tpn + sin ip n + d x d y p n = 0, 



r 



dlpn + Pn + d x dyip n = 0, 



(2) 



which reduces to the usual coupled sine-Gordon equa- 
tions for 7 2 3> 1. Note that a Lagrangian approach for 
SIJJs can be formulated only for two interacting fields ip 
and p, but not for tp alone. This because of the 2D nature 
of the vector potential in SIJJs. So particles with two 
types of polarization can propagate. For a ID Josephson 
junction, only one polarization is enough. 

Linearizing Eqs. J3J) results in the spectrum co 2 = 
1 + fe 2 /(l + ky) of the classical JPWs in the contin- 
uous limit (i.e., k y s <C 1) and 7 2 ^> 1. Here, k x 
and k y are the wave vectors (momentums in the quan- 
tum description; here, h = 1) of the JPWs. In or- 
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der to quantize the JPWs we introduce the Hamilto- 
nian, 7i — J2 n I dx(U Vn (p n + IL Pn p n ) — C, with the mo- 
menta and n Pn of the <p n and p n fields, and require 
the standard commutation relations \Q. Vn (x),ip. n i(x')\ = 
-iS(x - x')5 n ^, \n Pn (x),p n >(x')] = —i5(x - x')8 n ^ 
(all others commutators are zero), where 8 is either a 
delta function or Kronecker symbol. Expanding cos tp n — 
1 - ip 2 n /2 + ipi/24 - ... , we can write H = H a + H &a , 
where we include terms up to y> 2 in 7io, and the an- 
harmonic terms in 7i an . Diagonalizing TIq, we obtain 
the Hamiltonian for the Bosonic free fields a and b: 
= Ek v J( dk */ 27T ) {ea(fc) a+a + e b (k) 6+6)}, where 
the energy of the quasiparticles are 



£«(*) = 1 
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e b (k) = 
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up to I/7 2 , for 7 ^> 1. The energy e a (k) coincides 
with the frequency u(k x ,k y ) for classical JPWs, while 
the quantum bosonic field b corresponds to the gapless 
branch of the excitations. The original fields ip n , p n in 
Eq. are related to the free Bosonic fields a and b by 



zH+ + b 



r 



P ' 



b+ + b 



,(4) 



where Z = k x k y / (ky + 1). The case of a single Josephson 
junction corresponds to k y — resulting in a single field, 
o. 

Thermodynamics of quantum JPWs. — Finite temper- 
atures, T, excite both a and b quasiparticles provid- 
ing contributions of the JPWs to the thermodynam- 
ical quantities. The thermal equilibrium internal en- 
ergy E(T) = J2k v I dk xl ( 2?r ) [ £ a n a + £b rib] of the system 
can be calculated using the usual Bosonic distributions 
n a,b = l/[exp(e Qi b/7 1 ) — 1]. The calculated dependence of 
E(T) for SIJJs and, for comparison, for an "equivalent" 
stack of non-interacting Josephson junctions is shown in 
Fig. la. This clearly shows that the thermodynamic en- 
ergies are significantly different for these systems, espe- 
cially at low temperatures. Thus, finite temperatures 
easily thermally excite JPWs in layered superconductors, 
compared with the case of non- interacting junctions. The 
main origin of this enhancement is the suppression of 
both excitation energies e a ,b(k x , ky) when increasing k y , 
which is associated with a stronger interlayer interaction. 
Other thermodynamic quantities, e.g., heat capacity, can 
be easily calculated using the standard expressions. 

Interaction of Bosonic fields. — The interaction be- 
tween the a and b fields, including the self-interaction, 
occurs due to the anharmonic terms in 7i an w 
(— 1/24) J2 n Jdx fn +•■■, where <p is given in (@J. Here 
we consider the dominant first-order perturbation terms. 
Using the interaction representation, we obtain the am- 
plitude Sif = 2m(f\TC an \i)S(ei — £/) for the transition 
from the initial state \i) to a final state |/), where Eij 



are the energies of the initial and final states. To first- 
order approximation in I/7 2 , a decay (see Fig. lb) of an 
a-JPW can occur in two channels: either 3a or 2a + 6. 
The amplitude, Sdccay, °f the decay of the quantum a- 
JPW propagating along the £-axis, i.e., along the layers, 
k\ = (k x \, 0), is determined by 



fide 



,Y\d 2 ki \ s ( £ a(ki) -E 



2 8 3tt 5 



£a(kl) Y[ £ a(h) 



(5) 



5(e a (ki) - e a (k 2 ) - £ a (k 3 ) - £b(fei)J 



s(ki~Y.h 



(4 7 2 /9) Z 2 (k 4 )^Je a (k 1 )e a (k 2 )e a (k 3 )e b (k 4 ) . 

where the sums and products are performed over the final 
states with momenta ki. Eq. © predicts the probabil- 
ity I Sdccay 1 2 to create JPWs propagating perpendicular 
to the layers by an a-JPW quantum propagating along 
the layers. Using Eqs. Q and J5J, one can conclude that 
the amplitude Sdecay diverges for large k y i, when reso- 
nance conditions [e a (ki) — 2 or 3] are fulfilled. For the 
former case (in dimensional units, £ a (k\) — 2Tiujj) the 
a- JPWs create a-6-JPW pairs, while at e a — Shuij the 
2a-excitations diverge. Indeed, due to the ip 4 nonlinear 
interaction, a particle can only create two more addi- 
tional particles, which could be either 2a or a + b. The 
first process has a threshold 2fiujj (similar to the 2mc 2 
rest energy threshold for e~ + e + pair creation in QED), 
while the second one has a huij energy threshold due to 
the gapless nature of the b particles. Figure Id shows 
the calculated probability, | Sdccay | 2 , of decay of a JPW- 
a-quantum versus the energy e a of the initial a-quantum. 
Both resonance peaks are clearly seen. 

We can similarly analyze the scattering of a- JPWs. 
The diagrams in Fig. lc show two input particles as 
the initial state corresponding to particles "1" and 
"2", while the final state |/) contains free particles "3" 
and "4" . These diagrams do not diverge for any input 
particle momentum. However, the scattering probabil- 
ity enormously increases for large transverse momentum 
transfer (k y i — k v a), if the energies of the initial particles 
are close to fujj. This can occur either for low k x or 
large k y of the particles 1 and 2. The decay and scat- 
tering resonances occur due to the unusual anisotropic 
spectrum of the JPWs, i.e., e a (k Xl k y — > 00) = 1 and 



£b(k x , ky 



k x h- 



Enhancement of macroscopic quantum tunneling. — 
Now we apply our theory to interpret very recent ex- 
periments on MQT in Bi2212. To observe MQT, an 
external current J, close to the critical value J c , was 
applied 7J. This produces an additional contribution 
jip n in the Lagrangian Q). When tunneling occurs, the 
phase difference in a junction changes from to 2-7T, 
which can be interpreted as the tunnelling of a fluxon 
through the contact. This process can be safely de- 
scribed within a semiclassical approximation and we use 
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FIG. 1: (Color online) (a) Thermal energy (E)(T) of quantum 
JPWs versus temperature T for SIJJs (top red solid line) . The 
energy (E) is normalized by huij and sample volume. For 
comparison, the bottom dashed blue line shows (E)(T) for 
the same number of non-interacting Josephson junctions and 
with the same parameters. The dashed-dotted green and the 
dotted orange lines correspond to the contributions of the a 
and 6 fields to the energy, respectively. The chosen parameters 
are standard for Bi2212: s = 15 A, X ab = 2000 A, 7 = 600, 
and uij/2n = 150 GHz. (b,c) Feynmann diagrams for decay 
(b) and scattering (c) of JPWs quanta due to anharmonic in- 
teractions, (d) The probability per unit time, ISdccayl 2 , for an 
a- JPWs to decay, versus energy e a (ki) of the initial particle, 
for 7 = 300. Both axes are in dimensionless units. The blue 
dashed line in (d) corresponds to the a — > 3a channel pair- 
production and the red solid line to the a — > 2a + b channel. 



the approach developed in Refs. 0, to calculate the 
escape rate T = (tjp/27r)\/1207r.Bexp(— B) of a fluxon 
through the potential barrier. Here, top is the oscillation 
frequency of a fluxon near the effective potential mini- 
mum, and B = dr C(t = it) is described by the 
Lagrangian (JTJ with the classical fields determined by 
Eqs. (0), if we add the term j = J/J c in the right-hand- 
side of the first equation. In the limit 7 2 ^> 1, the equa- 
tion for (p is reduced to standard coupled sine-Gordon 
equations [||, which in the continuous limit, i; 9 s « 1 and 
y = ns/X a b, reads (l — d 2 /dy 2 ) [<p + sinp] — d 2 ip/dx 2 = 



j. We seek a solution of the last equation in the form 
ip = ip(x, y, t) + arcsin(j), where the field ip obeys 



1 - 



dy 2 



ip' — j(l — cosip) + \/l — j 2 sin-0 



d 2 ip 



= 0. 



3a; 2 

(6) 

Following the experimental setup pfl, here we consider 
the SIJJs having the size L ^> s along the y direction, 
i.e., the total number of contacts N = L/s 3> 1, and the 
size of the SIJJs in the x direction, 2d, is smaller than 
the Josephson length, Xj = 7-\/sA a &/2. 

We can linearize Eq. JHJ in all junctions except one, 
where the fluxon tunnels. The linearized equation can 
be solved by using the Fourier transformation, ip = 
E m / exp(-wjt) cos(k xm x)ip m (y,uj)duj/2ir, where k xm = 
A c 7r(2m + l)/2d. Since in the experiment [jj the sam- 
ple connects two bulk superconductors, we can choose 
the phase difference to be zero at the top (y — Li) and 
bottom (y = L\ — L) layers of the sample, and y = 
corresponds to the position of the fluxon tunneling. As a 
result, we derive the solution of the linearized equations 
in the form ip m (y) = ip m (G) smk[q m (Li - y)]/ sinh[g m Li], 
for y > 0, and ip m {y) = ip m (0) sivh[q m (L - L\ + 
y)]/ smh[q m (L 
,2 



Li)] for y < 0. Here, q 2 



(k 



2 

xm 



\J\ - j 2 - lu 2 )/(\/1 - j 2 - to 2 ). Following the method 
described in Ref. |2J and requiring the continuity of both 
ip and the current flowing through the central (y = 0) 
contact, we obtain the nonlinear equation for the phase 
difference in the junction with y = 0: 



dio 
2T 



\2 

ip — j(l — cos ip) + y/l — j 2 sinip = —j 



y-' k 2 xm sinh(g m ii) sinh(q m (L - Li)) 



sinh(g m L) 



cos(k xm x)ip m (u).(7) 



For an infinite (L, d — > 00) sample, this equation coin- 
cides with the nonlocal equation for the Josephson vor- 
tex in the SIJJs [2|. Eq. (Q can be used to describe the 
tunneling of a fluxon through a SIJJs with any width d 
and any number of layers N. However, such a treatment 
can only be done numerically. 

Now, we adopt Eq. J7J) for the short (d/Xj < 1) SIJJs 
used in 0, where d/Xj s» 2/xm / 5/im = 0.4. In this 
case, the phase difference tp changes slowly with x and 
the main contribution to the sum in the right-hand-side 
of JJJ) comes from the first harmonic k xm = irX c /2d. Ne- 
glecting contributions to the tunnelling process arising 
from higher- frequencies, to > — j 2 ) 1 ^ 4 , and integrat- 
ing Eq. (J2J) over dx, we derive for the phase difference 
ip, averaged over x, the equation: d 2 ip/dt 2 = —dV/dijj. 
Here the effective potential V(ip) can be written as 



ip 2 



V(i>)=j(smip-ip)-Vl-p(cosiP-l)-g n (j)^ , (8) 
where n = L\X a b/s labels the contact through which the 
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FIG. 2: (Color online) The MQT escape rate F versus di- 
mensionless external current i. Red and blue points are the 
experimental data, from Ref.0, for two different samples. Red 
and blue dashed lines are the curves To(j), Eq. 1111 . for these 
samples taken from 0. Red and blue solid lines are the func- 
tions r(j') calculated from Eq. Ill Oft , using data from Table 1 
of 0, for their samples US1 and US4, 7 = 600, s = 15 A. 
The inset shows the dependence of the escape rate V versus 
the number of contacts N in the SIJJs, using parameters for 
sample US4 of at j = 0.97 (red solid line). The blue dashed 
line shows NTq- 



fluxon tunnels, 



9n(j) = 



2(1 - j 2 ) 1/2 Q sinh(Qn) sinh[Q(iV 



sinh(QAT) 



(9) 



Figure 2 shows r(j), which very well describes experi- 
mental results in [7j. Some deviation between the ex- 
perimental data and the theoretical prediction at high 
currents is due to a significant lowering of the poten- 
tial barrier resulting in a decrease of the accuracy of the 
semiclassical approximation. The dependence of T on 
the number N of junctions is nonlinear due to the long- 
range interaction between different junctions, described 
by the last term in the expression (JHJ) for the effective 
potential. This nonlinearity is strong for relatively small 
N < N c — d/jL and the escape rate becomes propor- 
tional to N when the SIJJs thickness L exceeds the effec- 
tive interaction length d/j. Very different types of MQT 
models in SIJJs, with no quantitative comparison with 
experimental data, are also being studied in 12]. For 
instance, here we consider the inductive coupling among 
layers, which is known to be strong, instead of the weak 
capacitive coupling among layers used in |l2^. 

Conclusions. — We analyze the quantum effects in SI- 
JJs. We develop a model for quantum excitations in SI- 
JJs using two Bosonic fields. We also describe the inter- 
actions and thermodynamics of these fields. Moreover, 
we suggest a semiclassical theory of the fluxon quan- 
tum tunneling in SIJJs, which is in good agreement with 
recent remarkable experimental observations. The ob- 
tained results might be potentially useful for future de- 
signs of quantum THz devices. 
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and Q(j) — njs (l — j 2 ) ^ 4 /2d. For an applied cur- 
rent J close to J c , where tunnelling was observed 0, 
we can expand both cos^ and sin?/; and, finally, derive 
V$) = - Vi)/6, where = 3{^/)~p - 

9n(j)]- Using a semiclassical approach [lfj (i.e., B — 
Jq 1 [2V(t/>)] 1 / 2 d , ip) 1 and taking into account that the 
fluxon can tunnel through any junction of the SIJJs, we 
derive (now in dimensional units) 



r 
r 



N 



f- = >:(i 



\5/4 



n=Q 



exp 



36J7 
5 hup 



(10) 

where the summation is taken over all N contacts. Here, 
the effective Josephson frequency is u)p(j) = u)j(l — 
j 2 ) 1 / 4 , the height of the potential barrier Uq = 2Ej(l — 
j 2 ) 3 / 2 /3, the Josephson energy Ej = <& J c /2ttc, and the 
escape rate T (j) for a single Josephson junction (see, 
e-g-, 0) is given by 



Tot?) = 



TT V TlUJp(j) 



exp 



36 



Uo(j) \ 



5huip(j)J ' 



(11) 
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